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Criteria are developed for the practicability of expansion engines and expansion turbines 
in large hydrogen liquefiers. Some additional data on pertinent properties such as the sonic 
velocity of normal hydrogen are included. The analysis suggests that use of expansion 
turbines in hydrogen liquefiers having a production capacity of less than 5,000 liters per 
hour ordinarily would not be justified. For greater capacities, however, use of modern 
materials and design should make expansion turbine performance even more favorable 
than in liquefiers for denser fluids. 



1. Introduction 

New processes using large quantities of liquid 
hydrogen require the development of Liquefiers of 
much greater productive capacity than any existing 
plants. The magnitude of this increase introduces 
many engineering design problems and shifts em- 
phasis to parameters of lesser importance in lab- 
oratory production units. 

In particular, only thermodynamic cycles which 
result in a high overall plant efficiency can be con- 
sidered. Using the general properties of fluids and 
the history of the air liquefaction industry as guides, 
the course of efficient design can be charted with 
confidence; larger liquefier capacity will require 
increased use of expansion engines at certain stages 
of the cycle. 

The increasing interest in the use of such expansion 
devices suggests that the following analysis, originally 
issued as an informal communication in late 1957, 
might be of value to a larger audience. A few 
changes and additions have been introduced to 
make the ideas current. An excellent introduction 
to the use of expansion engines and turbines, in- 
cluding both historical commentary and descriptions 
of machines, has recently been published by Collins 
and Cannaday [l] 2 which also includes analyses 
of Land [2] and Jekat [3]. 

The type of expansion engine that should be used 
in a particular liquefaction plant depends, again, on 
the production capacity. Engines which have been 
developed sufficiently in other applications to be 
considered in hydrogen liquefier design include 
(listed roughly in order of increasing flow capacity) : 
reciprocating expanders; "positive displacement" 
rotary expanders, especially the Lysholm design; 
radial flow turbines; and axial flow turbines. 

Since liquefaction of hydrogen requires refrigera 
tion over a wide temperature range (20° to 300° K), 
expansion engines may not be the most desirable 
means of supplying such refrigeration at all stages 
of the process. At higher temperatures (200° to 
300° K), vapor-cycle units or vapor-cycle cascades 



1 This work was supported by the U.S. Atomic Energy Commission. 

2 Figures in brackets indicate the literature references at the end of this paper. 



using readily available refrigerants would be inex- 
pensive to procure and operate. Precooling to the 
temperature of liquid nitrogen (65° to 80° K, de- 
pending on the vapor pressure maintained) would be 
advantageous, since liquid nitrogen plants have been 
extensively developed, can readily be procured, and 
will operate at relatively low cost . 

At the lower end of the liquefier (due to the cooling 
of the gas to the temperature range 20° to 35° K) the 
yield of a simple Joule-Thomson expansion process is 
high. The overall plant yield a* ( = mass rate 
liquefied/mass rate through the compressor) can be 
only slightly increased by substituting an expansion 
engine. At this time, problems in the development 
of an expander, operating in a range 4 where the fluid 
condenses, more than offset such possible small gains 
in production. Thus, the expected operating range 
of expansion engines for hydrogen liquefaction lies 
within the region 30° to 80° K. Figure 1 shows the 
lower stages of a typical one-expander cycle. 

The low temperatures at which hydrogen ex- 
panders will operate may introduce new problems 
in the design of power-output devices. Clearly, the 
adiabatic extraction of power from a gas at some 
desired refrigeration temperature will be advan- 
tageous only if this power is transferred to some 
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higher — preferably ambient — temperature where it 
can be used (or dissipated, if it is not large) without 
being reintroduced to the gas as heat. 

Power-output problems are less important for low- 
speed engines since shafts or connecting rods can be 
long enough to make heat conduction negligible. 
Thus, such engines can be directly coupled with a 
variety of ambient temperature devices, including 
water-cooled brakes, oil pumps, fans, and motors. 
High-speed turbines, however, require short, thick 
shafts to maintain high critical speeds and thus 
would have higher heat flows if directly connected to 
ambient temperature output devices. More im- 
portant, high temperature-gradients would increase 
the danger of shaft warpage during operation. 

Since high-speed (75,000 to 250,000 rpm) turbines 
appear to be particularly applicable to hydrogen 
liquefiers in the 5,000 to 50,000 liter/hr range, 
practical solutions are desirable. Two immediate 
alternatives may be considered. The first, a direct- 
coupled high-speed hydrogen compressor operating 
in a liquid nitrogen bath, might be feasible if a high 
compressor efficiency could be achieved . The second , 
and thermodynamically more satisfactory method, a 
direct-coupled, two-pole, high-speed alternator driv- 
ing a multipole, low-speed motor at ambient temper- 
ature, is especially appealing, since electrical energy 
is "unaware" of temperature level, and since a simpler 
arrangement of liquefier components could be made. 
For such an alternator, the large decrease in con- 
ductor resistivity would result in much lower PR 
losses than at ambient temperature; the increased 
core permeability would permit higher flux densities 
and, thus, higher voltages; operation of the rotor in 
a low temperature hydrogen atmosphere (with low 
viscosity) would make rotor windage losses negligible. 
High efficiency would be important (since losses 
would be immediately returned to the hydrogen 
stream as heat), but the above temperature effects 
should allow exceptionally high performance. 

2. Types of Expansion Engines 

Reciprocating expansion engines have been used 
for years in both liquefiers and low-temperature 
process plants. Crankshaft-type engines (in con- 
trast to walking-beam) are used in the Amarillo 
plant of the Bureau of Mines for the quantity separa- 
tion of helium from natural gas. Engines with very 
long, thin connecting rods operate to 10° K in the 
commercially produced Collins' cryostat, and in the 
45 liter/hr helium liquefier constructed by Collins 
for the Massachusetts Institute of Technology 
Cryogenic Engineering Laboratory. Application of 
reciprocating expanders in medium capacity hydro- 
gen liquefiers should involve only size and details 
of design. 

The upper limiting liquefier capacity for which 
reciprocating engines can profitably be used depends 
on the bulk which the plant designer is willing to 
insulate and maintain at the required temperature 
level. Therefore, this limit must depend to some 
extent on good engineering judgement. For the 



above reasons, this study was restricted to problems 
in the application of rotary expanders. However, a 
brief consideration of such possible upper limits for 
reciprocating expanders indicated that their use 
remains practical in liquefiers up to a capacity of at 
least 10,000 liters/hr. 

For other reasons — lack of data — the possible 
application of the Lysholm "positive-displacement" 
machine has only been noted. A compressor built 
by the Svenska Rotor Maskiner AB has been tested 
as an expander for air at pressure ratios from 3 to 
9 and flows from 8 to 28 lb/min. A graph showing- 
resultant test efficiencies (on a temperature, rather 
than enthalpy basis), furnished by the Stratos 
Division, Fairchild Engine & Airplane Corp., is 
included in figure 2 of this report. Maximum 
efficiency does not appear to be exceptionally high — 
approximately 73 percent for pressure ratios 6 to 
9 and flows from 20 to 28 lb/min at a speed of 40,000 
rpm. Since a large part of the losses in such machines 
can be ascribed to leakage, their application to hy- 
drogen liquefiers is uncertain. 

Turbine expanders may be classified by flow di- 
rection and type of blading. Three main flow types 
are: axial flow, radial inward flow (centripetal tur- 
bine), and radial outward flow (centrifugal turbine). 
Radial outward flow turbines have smaller pressure 
and enthalpy drops and somewhat poorer efficiencies 
than radial-inflow turbines. Probably, an outward 
flow turbine would only be advantageous in a plant 
without a final J-T stage, since cycling of fluid 
that might condense within a stage would be pre- 
vented. 

The choice between axial-flow and radial-inflow 
turbines depends on a number of complicated factors, 
including required volume flow and pressure ratio 
and allowable efficiency. Because of the centrifugal 
field, a radial stage will permit higher pressure and 
enthalpy drops, and may be as much as 10 percent 
higher in stage efficiency (78 percent is usually con- 
sidered a maximum for axial stages: efficiencies 
of 90 percent are not uncommon for well-designed 
radial stages). Pressure ratios as high as 17 have 
recently been obtained in a commercially produced 
supersonic radial inflow turbine with total peak 
test efficiency remaining over 72 percent (supersonic 
flow will ordinarily result in lower efficiency). 

Four major losses exist in a turbine stage: nozzle 
flow loss, rotor flow loss, leakage, and disk friction. 
In small mass-flow machines operating on large 
heads, leakage may be a decisive factor in deter- 
mining the efficiency, especially with hydrogen. The 
leakage loss in an axial-flow stage may be as high as 
nine times the leakage loss of a radial stage. 

If the required pressure ratio is such that a few 
radial stages can be used, better performance can be 
expected. For very high pressure ratios, however, 
despite the smaller stage drop, an axial-flow turbine 
would be preferable since there would be no turning 
losses such as are associated with the repeated change 
from axial to radial flow in multi-stage radial tur- 
bines. Configurations of typical axial-flow and 
radial-inflow stages are shown in figure 3. 
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Figure 2. Performance chart for Lysholm expander. 





AXIAL-FLOW STAGE RADIAL INFLOW STAGE 

Figure 3. Typical r<t<li<tl ami axial flow stages. 

Dimensional analysis can be used to give some in- 
dication of the proper choice of flow-type for a partic- 
ular application. Such an analysis shows that an 
important dimensionless parameter in turbine design 
is the specific speed 7i s =n-\W/p(AH°) 5/2 , where n is 
the actual speed in rps, W the power output in watts, 
p the gas density in kilograms per cubic meter, and 
AH the isentropic enthalpy change in joules per 
kilogram. Using the relation W=7]firhAH°, the ex- 
pression can be transformed to n s =n-\ vfiV/ (AH°) 3/2 , 
where the actual volume flow /3V through the turbine 
is usually specified at outlet conditions. Thus, with 
equal speed and head, larger volume flows imply 
greater specific speeds. 

The above equation for specific speed can be ex- 
pressed in terms of a number of velocity and diameter 
ratios. By varying these ratios in a range about 
their optimum values, the performance of any type of 
turbine can be determined graphically as a function 
of specific speed if reasonable values of certain loss 
factors are assumed. Such a calculation has been 
made by von der Nuell [4] and leads to the con- 
clusion that, up to a specific speed of about 0.1, the 
radial-inflow turbine is preferable to the axial-flow 



turbine. One of his graphs, showing ranges for n s as 
a function of the velocity ratio r v (the ratio of the tip 
speed u and the isentropic head velocity c =^2AH°) 
is given here in adapted form as figure 4. (Note that 
in consistent units, n s should be dimensionless ; how- 
ever, von der Nuell's equations are in the English 
gravitational system using hp, rpm, and head in feet, 
while his graphs are in a metric gravitational system 
using metric hp, rpm, and head in kilogram-meters 
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Figure 4. Ranges of application for 50 percent reaction 
axial-flow turbines - - -, impulse (R — 0%) turbines 

— . . — . . — . — ^ anc i fig percent reaction radial-inflow 
turbines — ■ — ■ — , in terms of specific speed n s and velocity 
ratio r v = u/c for practical diameter ratios 5 = d/D and 
nozzle angles a. (From von der Nuell.) 
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per kilogram; thus his equations give values of n s 
which are 0.9 times the values shown in his graphs; 
in this paper, all values are in the absolute mks sys- 
tem and, thus, will be the same as values in any other 
consistent absolute system, but will be 0.026 times 
the values in his graphs.) 

As an example of the specific speed concept, cycle 
calculations [5] indicate that, for expander inlet con- 
ditions Pi=30 atm, 7^ = 65° K, an optimum mass 
flow ratio is 0=0.46. With an expander efficiency 
of only 80 percent an overall yield of a* = 0.35 should 
be obtained. The theoretical head would be AH°= 
429,000 j/kg, the isen tropic head velocity c =927 
m/sec, and, at actual outlet conditions, the density 
would be 13 Amagats. In terms of liquefaction ca- 
pacity, C Ll in liters per hour, the specific speed for 
these conditions can then be expressed as n s =0.5X 
10- 6 <nJCT m 

If n s = 0.l is taken as an approximate upper limit 
for the practicable application of radial inflow ex- 
panders, then the limiting liquefaction capacity for 
this case is 6\ = 4X I0 10 /n 2 . The highest operating 
speed for which standard antifriction bearings are 
commercially available is approximately n = 2,000 
rps (120,000 rpm), which would permit a maximum 
liquefaction rate of 10,000 liters/hr. Similarly, for 
a more conservative speed of 60,000 rpm, the maxi- 
mum rate would be 40,000 liters/hr. 

A similar calculation can be made using the same 
conditions, but a minimum specific speed of approxi- 
mately n s = 0.025. Here, it should be emphasized 
that many other factors must be considered in a 
detailed design, such as minimum practicable 
dimensions and the increasing percent tolerances 
that must be allowed with decreasing dimensions. 
However, using specific speed merely to obtain an 
approximate range of applicability, the lower limit 
for liquefaction capacity can be expressed as C L{m]n) = 
25X10 8 /ft 2 . For a 120,000 rpm turbine, this mini- 
mum capacity would be 600 liters/hr; and, for 
a 60,000 rpm turbine, 2,500 liters/hr. Minimum 
liquefaction rates are estimated more exactly in 
section 7, and in addition, the effect of the use of 
cold gas bearings (permitting speeds up to approxi- 
mately 250,000 rpm) on the minimum liquefaction 
rate is evaluated. 

In conclusion, the entire range of foreseeable re- 
quired liquefaction capacity can evidently be covered 
by two types of expansion engines: reciprocating 
expanders for small to medium capacity, and radial- 
flow turbines for medium to large capacity. The 
radial inflow turbine will have the advantage of high 
stage enthalpy and pressure drops, and high effi- 
ciency. A further advantage is that a given turbine 
can be adjusted for a large range of flows, without 
sacrificing efficiency, by the use of adjustable nozzle 
blades — a frequent practice in small radial-flow air- 
craft turbine design. Such variable flow character- 
istics are desirable to cover uncertainties in the 
original cycle design and fluctuating conditions 
during operation. 

The application of radial turbines requires solu- 
tions to several problems: (a) Proper design of the 



turbine, itself, for the conditions imposed by the 
cycle; (b) determination of the important fluid 
properties, and, for the particular case of hydrogen, 
their effect on output and performance; (c) proper 
design of a power-output device which, since the 
radial turbine is a characteristically high-speed 
machine, might take the form of a high-speed 
alternator; and (d) proper design of turbine (and 
generator) bearings which — again, because of the 
high speeds — might be required to operate at low 
temperature. 

3. Fluid Properties 

Values of the density and enthalpy of hydrogen 
required for turbine design can be obtained from the 
paper of Woolley, Scott, and Brickwedde [6]. In 
order to determine the compatibility of the turbine 
expansion stage with the J-T stage, curves were 
calculated (and extrapolated at the extremes of 
pressure) for the yield a as a function of pressure 
with the inlet temperature to the J-T heat exchanger 
as a parameter. These curves, determined for heat 
exchanger efficiencies of 100, 95, and 90 percent using 
the method described by Brown and Dean [7], are 
given in figures 5, 6, and 7. For an efficiency of 100 
percent, the low temperature limiting curve was 
obtained by integrating the expression (dh/dp) T = 
V-T(bV/dT) p using the liquid PVT data of Fried- 
man [8]. Higher isotherms were obtained from the 
enthalpy data of [6] up to about 120 atm. The re- 
mainder is based on extrapolation and general 




Figure 5. Yield a of the J-T stage for a heat exchanger 
efficiency, i)=100 percent for normal hydrogen. 
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Figure 7. Yield a of the J-T stage for a heat exchanger 
efficiency, v = 90 percent for normal hydrogen. 

thermodynamic tendencies and, thus, should only be 
considered to show the general behavior of the yield 
curves. These curves emphasize an important 
design factor: for small liquefiers using turbine ex- 
panders it is desirable to keep the expander inlet 
pressure low enough thai reasonable flow areas can 
be maintained; the opt i mum pressure for the J-T 
stage alone will ordinarily be much higher; 1 1ms, some 
compromise is required, preferably a pressure on, or 
just above the "knee" of the yield curve. If the 4 
temperature is low enough, however, such pressures 
should require only two, or, at the most three, radial 
stages in a subsonic turbine, and one or two stages 
in a supersonic turbine. 

In addition to density and enthalpy, the sonic 
velocity c and kinematic viscosity v are properties 
of primary importance in turbine design. The 
stage work that can be produced under comparable 
conditions is directly dependent on the sonic velocity. 
Similarly, the major part of the stage losses— and, 
thus, the efficiency — is a function of the kinematic 
viscosity. 

Only a small amount of data is available on the 
sonic velocity of hydrogen. NBS Circular 564 [9] 
gives values at 0, 1, 10, and 100 atm, but includes 
only a few points below 80° K. For this reason, a 
detailed calculation of the thermostatic sonic ve- 
locity c=l/V pk s was made from the identity 

Mc 2 /d P Z\ ,R/dTZ\ 2 
RT V ^ J T ^c v V dTj p 

using the PVT and specific heat data of [6], where R 
is the gas constant, M the molecular weight, T the 
absolute temperature, p the density, c v the specific 
heat at constant volume, and Z the compressibility 
factor. The results are given by the combination 
of the results in table 1 of the sonic velocity in the 
perfect gas state c° = ^RTy°(T)/M and the ratio 
c(p, T)/c°(T) given in figure 8 as a function of tem- 
perature at densities from to 500 Amagats (units of 
density reduced by the density at standard conditions 
of 1 atm and 0° C). 

The isen tropic critical velocity and critical pres- 
sure ratio of a nozzle can be obtained from the 
curves of figure 8 together with a T-S chart by an 
approximation process. Maintaining a constant 
entropy corresponding to the initial conditions P u 



Table 1. Sonic velocity c°(T)=J~ Ty°(T) of normal 
hydrogen in the perfect gas state 
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Tf, the final value of the static enthalpy //,- can be 
estimated. This estimate would result in a final 
velocity v f =-\ l 2(h i —h f ) which can be compared with 
the sonic velocity from figure cS based on the same 
assumed P f and T f . By adjusting the final state 
until v f =c(P f , T f )=v c , both the critical velocity and 
pressure ratio are determined. 

Similar calculations of the kinematic viscosity v of 
hydrogen were made using data for the absolute 
shear viscosity and density from [6,8,9], and some 
unpublished data of D. B. Chelton et al., of the 
National Bureau of Standards, Cryogenic Engineer- 
ing Laboratory. The results, in figure <), give the 
kinematic viscosity v in square meters per second in 
the range 14° to 100° K at pressures up to 100 atm. 

4. Comparison of Different Fluids 

In order to determine the effect of differing fluids 
on turbine output and performance, some reasonable 
basis for comparison must be assumed. Since the 
final goal is liquefaction, it appears reasonable to 
compare only turbines at similar points in similar 
cycles in plants having the same liquefaction ca- 
pacity in kilograms per second. In addition, it is 
clear that the comparison of a hydrogen expander 
and a nitrogen expander at the same inlet tempera- 
ture would not be significant: 7^ = 60° K would be 
reasonable for hydrogen, but at this same point 
nitrogen would be a solid. To obtain a practical 
basis for comparison, reduced quantities can be 
defined in a number of ways, one of which is 
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Pc—Pb 



where c and b refer to the critical point and normal 
boiling- point. This set is particularly suitable in 
that it relates the temperature and pressure to the 
values at which liquid is to be produced as well as 
to the critical point values characteristic of the 
particular gas. 
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Figure 8. *Sonfc velocity of normal hydrogen reduced by perfect gas values. 
Note: A large size copy of figure 8 is available upon request to Cryogenic Data Center, National Bureau of Standards, Boulder Laboratories Boulder, ( 
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Figure 9. Kinematic viscosity of normal hydrogen. 

The purpose of the expander is to supply refrig- 
eration to a heat exchanger so that the gas to be 
liquefied can be cooled sufficiently for introduction 
to the J-T stage. However, temperature drop is 
not significant by itself, but must be related to a 
temperature level: e.g., a 1° K drop in a gas initially 
at 1° K would be impossibly large. Use of the 
reduced temperature r provides a meaningful meas- 
ure of cooling. For small increments, equal values 
of At = t AT/T can be considered to signify equiva- 
lent cooling for different gases. 

Using this measure, enthalpy balance for the first 
and second streams of a heat exchanger for a par- 
ticular gas a requires that 

where mean specific heats and temperature levels 
over each stream are used as an approximation 



which is sufficiently justified for the present purpose- 
determining a practical basis for comparison of 
turbines using different fluids. To the same degree 
of approximation, the mean temperature levels of 
each stream can be considered the same, T\~T\. 
Writing a similar equation for another gas b, and 
combining the two equations gives 



m\lm\ c p \/c» 2 At?/At» 

For equivalent cooling, At\=At\ and At^AtJ, and 
since the specific heat curves of all gases have a 
common characteristic behavior, at the same re- 
duced pressure and temperature ~c a v x[c a n will be 
approximately equal to cji/cj 2 - Thus, the follow- 
ing equation, 

t 6 






m\ 



must hold. Stated simply, the fact that one re- 
frigerating gas has a higher specific heat than anol her 
does not change the required mass flow ratio, since 
a gas of correspondingly higher heat capacity must 
be cooled. 

To make the transition from refrigeration mass 
flow to turbine mass flow, a second assumption 
that for the same reduced pressure and temperature 
the yield a of the J-T stage will be approximately 
equal for different gases — must be introduced. A 
comparison of figures 5, 6, and 7 with the similar 
curves for helium in [7] shows that this assumption 
is good: for example, the limiting pressure for hy- 
drogen is roughly 400 atm. the limiting pressure for 
helium is 42 atm; thus, limiting reduced pressures 
are tt=(42-1)/(2.3-1)=31.5 for helium, and w= 
(400-1) /(13-1) =33 for hydrogen. Together with 
the previous equation for refrigeration mass flow, 
this assumption requires that & n rh a =fi l) 7h b ] that is, 
for similar liquefiers with the same mass rate of 
production of different liquids, turbines at the same 
reduced pressure and temperature levels must have 
the same mass flow as each other to give equivalent 
refrigeration. A practical basis for comparison of 
similar turbines using different gases has thus been 
obtained. 

5. Turbine Output 

The theoretical output of a turbine stage can^be 
obtained from the Euler turbine equation! AH= 
A(u-v), or 

H i — H f =u 1 v l cos a\ — U2V2 cos a 2 , 

where 1 and 2 refer to conditions at the inlet and 
exit of the rotor, u is the tip speed, and v the ab- 
solute gas velocity. The inlet gas velocity v can 
be assumed approximately equal to the nozzle exit 
velocity — that is, having the same magnitude and 
direction, with the direction, also, approximately 
equal to the actual nozzle blade angle. To obtain 
maximum work, the rotor blading is ordinarily 
designed to make the exit whirl as small as possible. 
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With such a design, the second velocity term in 
the above equation can usually be neglected. Thus, 
the theoretical output is given approximately by 

\AH\=U V Y X . 

Since the Euler turbine equation in the above 
form is already in terms of mass flow (the units of 
AH could be watts per kilogram per second), one 
of the requirements for comparison of turbines with 
different fluids is already satisfied. Therefore, 
similar turbines having identical tip speeds and inlet 
velocities must have the same theoretical output for 
all fluids. 

The assumption of identical nozzle exit velocities, 
however, cannot be made if similarity in the flow is 
to be maintained. A value of i\ which indicated 
subsonic nozzle flow for hydrogen might require 
supersonic flow for nitrogen. To avoid this difficulty, 
the Euler equation can be transformed by defining 
a nozzle Mach vector Mi^v^c, where c is the sonic 
velocity for the actual static temperature and pres- 
sure which occur at the throat, and an overall 
"turbine Mach vector" M* = Ui/c, giving 

|Aff|=<MK*.Mi 

If the vectors M* and M^ are maintained constant 
with different gases the flow conditions and velocity 
triangles will remain the same. For sonic conditions, 
where Mi=|Mi|=l, since (dG/dp) s =0, the mass 
flow per unit area G will be maximum and c will be 
equal to the nozzle critical velocity v c . Since 
(dc/dp) s >0 the local sonic velocity c will always be 
less than the sonic velocity c f for the initial conditions 
P iy T t . However, c and v c will always be increasing 
functions of c<, so that gases with higher values of 
c t will also have higher values of c and v c . In terms 
of reduced pressure and temperature the Euler 
equation can then be expressed as 

where <j> is a function too complicated to be evaluated 
other than by approximate methods; however, for 
the approximation needed for this comparison such 
evaluation is not necessary. 

The last equation shows that more power output 
can be obtained under similar conditions from a 
turbine operating on hydrogen than any other gas, 
since the sonic velocity of hydrogen at the same re- 
duced temperature and pressure is greater than for 
any other gas. As an example, for the conditions 
in* section 2, T* = 65° K, P t =S0 atm, the sonic 
velocity of hydrogen is 680 m/sec. At similar 
conditions (r*=2.38, 7r=2.42), the inlet pressure 
and temperature for a nitrogen turbine would be 
P t = 79 atm, T*=249° K, and the initial sonic 
velocity, c* = 336 m/sec. Thus, the hydrogen turbine 
would have approximately four times the output 
of the nitrogen turbine. 



6. Turbine Losses 

The four main losses in a well-designed turbine 
operating near its design point are: Nozzle flow 
losses, rotor flow losses, disk friction, and leakage. 
For a turbine operating awa}^ from its design point, 
secondary flow losses may become as great or 
greater than all other losses combined. In addition, 
a supersonic turbine may have large shock losses 
which would be very difficult to predetermine 
analytically. In this discussion, however, only losses 
at the design point will be considered. 

An exact determination of nozzle flow losses would 
include effects due to the compressibility of the fluid 
and the change in flow area. Similarly, rotor flow 
losses should also include effects of turning and the 
centrifugal and Coriolis fields. Despite the probable 
importance of these effects, however, it is common 
practice in turbine analysis to avoid the formidable 
difficulties presented by them and base estimates of 
turbine losses on semiempirical equations such as 
the Prandtl-Karman law for the flow of incompress- 
ible fluids in stationary straight pipes of constant 
cross section. Such estimates have been found to be 
in approximate agreement with the results of actual 
tests [10] and also [1, 4j. The drop in total pressure 
P over a mean flow path length L is given by 

when D is a mean hydraulic diameter. This equa- 
tion can be considered, not as a physical law, but 
merely a definition of the dimensionless resistance 
coefficient X. Since such a loss in total pressure 
represents a conversion of mechanical energy into 
unusable heat, the total enthalpy loss XH can be 
determined from 8H=AP/p, that is 



8H=\-^-\v>=\-^-\ Ci ?-M> 



(0 



Defining \p=Co/c t and using <j> = c/c t as before, the 
relative loss can then be written 



SH , L ,,, ^ v 



(IT 



where the loss coefficient X is given by the Prandtl- 
Karman law [11] 

-^==2 log (Re A /X) -0.80. 

-y X 

Graphical solutions of the above equation show that 
X is a monotonically decreasing function of the 
Reynolds number. Since He=vD/v=MD c/v can be 
expressed in the form 



lie 



-v 



mMc 

irpv' 



the assumption of equal mass flow and Mach number 
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for different gases requires that the loss coefficient 
decrease as the factor \cjpv 2 increases. Comparing 
hydrogen with other fluids, the kinematic viscosity v 
at similar reduced pressures and temperatures will be 
approximately constant — at least for diatomic gases, 
while the sonic velocity c will be greater and the 
density p less. Thus, the flow losses in a com- 
parable hydrogen turbine can be expected to be 
somewhat less. 

Disk friction losses can be treated in a similar 
way. The moment acting on both sides of a disk 
of radius R rotating with an angular velocity co is 
given bv [12] 

1 

2 



^/=\*.^pco 2 .£ 5 



Using the specific speed at 100 percent efficienc}^ 
n° s , the relative disk friction loss can be expressed as 






M* 



(t)\ 



r%/2 n? 



The value of the rotational loss coefficient X* 
depends on whether the turbine Reynolds number 
Re* = / u 1 D 1 /2v = cM*D l /2v is in the laminar or turbu- 
lent range. For Re*<2000, X* is given by 



X*= 



h # Re* 



where h is the spacing between rotor and casing. 
For large Re*, the value of X* is independent of the 
spacing, and can be conservatively estimated from 
the Goldstein law [13] for a free disk 

-L=1.97 1og (Re*Vx*)+0.03, 
yX* 

which is identical, except as to constants, with (he 
Prandtl-Karman law. 

Although in either case X* is a monotonically 
decreasing function of Re*, a comparison of Re* for 
different gases cannot be made directly, since the 
rotor diameter D l cannot be immediately replaced 
by an expression involving mass flow. However, the 
larger volume flow in a comparable hydrogen turbine 
will require a larger rotor diameter to maintain good 
rotor passage design. Since c/v will also be greater, 
the disk friction losses can be expected to be less. 

The volume flow leakage of hydrogen will be much 
greater than other gases, since for comparable con- 
ditions the absolute shear viscosity will be much 
smaller: for r=2.38, w=0, the viscosity of hydrogen 
is /x=3X10~ 6 dekapoise (= kg /meter-sec) while the 
viscosity of nitrogen is n= 15 X 10~ 6 dekapoise. How- 
ever, it is the loss of mass from the primary stream 
rather than volume which must be used to determine 
the leakage loss. The pressure drop through a leak- 
age orifice can be written 



where the pressure coefficient X, assuming laminar 
flow, is given by X=64/Re. Solving this equation 
for the leakage mass flow m L in terms of reduced 
pressure drop gives 



m L 



Aw 1 



4irD l 7T V 



where A L is the leakage area, and \lTro = P c —P b . 
Both I/to and 1/v will be smaller for hydrogen than 
other gases (except helium) which will tend to offset 
the unfavorable increase in A L 2 /D X with the larger 
hydrogen turbine; thus, with careful design of seals 
and leakage passages it should be possible to main- 
tain equivalent leakage losses. 

A general conclusion can be drawn from the above 
discussions. The formulas used are somewhat over- 
simplified: an actual turbine designer would, for air, 
use more complex equations based on test data, with 
appropriate corrections in empirical constants. How- 
ever, if the design analysis used for air turbines, 
together with all the factors based on air- tests, are 
used for a comparable hydrogen turbine, the results 
should be reasonable and the performance actually 
better than indicated. 

7. Turbine Design 

After the above considerations it is possible to 
sketch a tentative design for a turbine in a specific 
application. Using as an example the conditions 
stated in section 2, that is, T,= 65°K, P,=30 atm, 
^=138 Amagats, c*=680 m/sec, and A,= l, 100,000 
j /kg, an expected efficiency of 80 percent would result 
in the outlet conditions 7 7 / =23.5°K, P/=l atm, 
P/=13 Amagats, and ^,= 756,000 j/kg. The theo- 
retical required isentropic head is A£P = 429,000 j/kg 
and the isentropic head velocity is Cq=^2AH° = 927 m/ 
sec. If the nozzles were designed for M=l, the criti- 
cal nozzle velocity would be v c = 4:80 m/sec and the 
critical pressure ratio would be X c =3.53. An addi- 
tional velocity is of some interest : the maximum the- 
oretically possible nozzle velocity, which in this case 
is v m =*\'2h i = 1483 m/sec, which can occur only for 
an infinite Mach number and pressure ratio. 

The first question that should be answered is 
whether the required conditions can be met with a 
single stage turbine. This can be done using the 
Elder equation assuming negligible outlet whirl, 
AH°=cl/2=UiViQO&oii. As a first possibility, as- 
sume a tip speed 1^ = 500 m/sec, which would not 
be considered unreasonable in the small aircraft tur- 
bine industry, and a minimum nozzle angle of 7° 
(for an axial-flow turbine a larger minimum nozzle 
angle would have to be assumed, since in that case 
the nozzle angle is with respect to a plane surface 
rather than a cylinder). These conditions result in 
an outlet nozzle velocity z>=875 m/sec, and a nozzle 
Mach number of M = 2.5 which must be considered 
to be much too high if the desired efficiency is to be 
obtained. Thus, using conventional rotor materials 
and limiting tip speeds, the desired head cannot 
reasonably be obtained in a single stage. 
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The possibility of newly developed materials, how- 
ever, tempers this negative answer. Recent tests on 
specially prepared aluminum alloys indicate tip 
speeds as high as 915 m/sec (=3000 fps) may be 
obtained in medium sized rotors [14]. Using the 
same nozzle angle as above and assuming 90° blades 
for greatest strength (and simplest fabrication), the 
equation c\=2uiVi cos a x becomes cl=2vl cos 2 a x and 
results in a tip speed of ^i = 655 m/sec (which 
should be reasonable with these newer materials) 
and a nozzle outlet velocity of ^=660 m/sec. At 
the nozzle outlet, the Mach number would be M= 1.3, 
which would not result in unreasonable losses. 

Since it is at least possible to design a single-stage 
turbine for the above conditions, it is of considerable 
interest to determine more precisely the minimum 
capacity liquefier which might use such a turbine. 
Again neglecting outlet whirl, the Euler equation can 
be written 



AH° 



nfim cot a x 



6i 



Pi 



where &i is the blade depth at the tip. To make the 
results most favorable, extreme conditions of 71=2,000 
rps (commercially available in an 8-mm ball bearing) 
and a probably minimum machinable value of 
hi — 2 mm will be chosen. To obtain the required 
tip speed of Ui = 655 m/sec, the rotor diameter would 
be approximately 10 cm. With the calculated nozzle 
velocity ^ = 660 m/sec, the density would be 75 
Amagats, but the assumption of a 3-percent relative 
nozzle loss would reduce this to an actual value of 
65 Amagats. With these figures, the minimum total 
volume flow at standard conditions becomes 7.38 m 3 / 
sec (=15,600 scfm), corresponding to a minimum 
liquefaction capacity of approximately 9,000 liters/1 
hr. 

Gas bearings have recently been the subject of in- 
tensive development. Air bearings for service up to 
250,000 rpm have been developed for commercial 
high speed grinding machines, and it seems clear 
that such bearings should be adaptable to low tem- 
perature service with hydrogen gas. Considering 
the use of such gas bearings, and choosing ?i=4,000 
rps, with the other thermal conditions as in the 
above example, would result, for a 5-cm rotor, in a 
minimum liquefaction rate of 4,500 liters/hr. Since 
a number of parameters have been chosen at ex- 
treme limits to obtain this result, it would appear 
more practicable to use a turbine of more than one 
stage or a reciprocating expander for liquefiers of 
smaller capacity. 

8. Bearing Tests 

The combination of high tip speeds required to 
make use of the enthalpy head of hydrogen, and 
minimum practicable flow areas for turbines used in 
small liquefiers, make the necessity of relatively high 
rotational speeds evident. Such speeds introduce a 
number of problems in bearing design, especially 
since, for several reasons, it is desirable to keep the 



bearings cold. Tanza [15] has found that small stain- 
less steel ball bearings with Micarta retainers will 
operate satisfactorily at 10,000 rpm without lubri- 
cation when surrounded by a cooling and inert 
(hydrogen) atmosphere. 

In order to determine whether similar bearings 
would operate at higher speeds, a small turbine-driven 
bearing tester was constructed. This was designed 
to test 10-mm bearings at 100,000 rpm at low tem- 
peratures. The turbine was an impulse type using a 
rotor that could be easily machined, and mounted 
on the same shaft as the pair of bearings being 
tested. The rotor was 1 in. in diameter, made of 
24ST-4 aluminum, and pressed on a hard-chromed 
Inconel shaft. Bearing races and balls were 440C 
stainless steel, held in a 440A stainless steel cylindri- 
cal block containing an adjustable spring to allow 
various axial loadings. The entire mechanical as- 
sembly fitted into a brass case of which the nozzle 
ring was an integral part. This case was designed 
to mount on the end of a Collins' ribbon-packed heat 
exchanger, and at the lower end had a filter to remove 
water or oxygen particles from a liquid nitrogen bath. 
The bearings themselves operated in nitrogen gas, 
but were cooled (and possibly partly lubricated) by 
liquid nitrogen jets directed between the races. A 
cross section of the tester is shown in figure 10 and 
several views of the tester and components in figures 
11, 12, and 13. 

Tests were completed on only one type, a bearing 
manufactured by the Barden Corp. using composi- 
tion retainers. Preliminary runs at room tempera- 
ture with oil-mist lubrication were highly successful, 
top speeds of 130,000 rpm being attained, and resulted 
in only one minor modification: the addition of 
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Figure 10. 



Cross section of high-speed low-temperature bearing 
tester. 
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Figure 11. Bearing tester: outer case 











Figure 12. Bearing tester: mechanical assembly. 




Figure 13. Bearing tester: assembly of rotor, shaft, and 
bearings. 

stiffening rings to the outer brass ease to raise its 
fundamental four-nodal frequency above 100,000 
rpm — exactly the designed operating speed. 

At low temperature (liquid nitrogen) the bearing 
tester itself operated just as satisfactorily, but the 
bearing retainers (which had appeared very promis- 
ing from Tanza's similar low speed tests) failed in 
apparently brittle fracture at the ball holes. It is 
assumed that the different results are due to the 
greater impact at high speeds of the balls oscillating 
about their mean rotating positions. The possibility 
of a difference in the low temperature brittle be- 
havior of the retainer material should be investigated 
furt tier, however, by standard notch impact methods. 
The highest speed attained at low temperature was 
60,000 rpm with one retainer cracked completely 
across at one hole. At this speed centrifugal force 
was great enough to deflect the broken retainer and 
stop the turbine until the retainer could return to its 
normal shape. In most tests retainer damage was 
considerably greater, large fragments frequently 
being broken completely away. 

Three conclusions can be drawn from the tests: 

(a) A satisfactory bearing tester for operation at 
100,000 rpm and 80° K or lower has been constructed; 

(b) bearings with composition retainers are probably 
not suitable for the combination of low temperatures 
and high speeds; and (c) microscopic examination 
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of the ball tracks, which were straight and extremely 
faint, indicates that ball bearings can probably be 
operated satisfactorily without lubrication at high 
speeds and low temperatures if supplied with inert 
coolants, provided a suitable retainer material can 
be found. In regard to the last conclusion, the 
manufacturer of the bearings tested has stated that 
similar bearings can be supplied with silver retainers. 
For this application silver would have the advantage 
of high ductility, low friction, and (unlike stainless 
steel) high thermal conductivity, and, thus, would 
appear quite promising in any continuation of such 
tests. 



9. Generator Design 

From the cryogenic point of view, the design of 
a high speed alternator as a possible power-output 
component for a turbine expander presents prac- 
tically no problems; conductor resistivity and core 
permeability are known sufficiently accurately for 
any temperature at which the generator might 
operate that the designer can merely obtain the 
proper values and thus eliminate temperature level 
from further consideration. Thus, most of the major 
problems, such as high speed, are not new and have 
been at least partially solved for room-temperature 
operation; a 90,000 rpm permanent field alternator 
is in commercial production and an 120,000 rpm 
induction motor produced by the General Electric 
Co. could, according to the manufacturer, be easily 
modified for use as an induction generator. 

Unlike present applications, however, the use of 
such generators with turbine expanders at low 
temperatures would require high generator efficien- 
cies since losses would be returned to the system 
as heat. With one exception, all such losses can 
be calculated with existing data and will apparently 
be much lower than at ambient temperature. 

The exception--hysteresis loss in the core — has 
evidently been of little interest: tests by Fleming 
and Dewar [16, 17] in 1896 showed no change in 
hysteresis loss in transformer iron at liquid air 
temperatures; tests by Terry [18] in 1910 indicated 
a rapidly increasing loss in electrolytic iron with 
decreasing temperature, but were not made at con- 
stant flux densities and thus are of questionable 
significance. 

Due to the scarcity and conflicting information on 
hysteresis loss, a modified "Epstein-test" was per- 
formed on silicon steel transformer laminations. The 
detailed results are reported separately by Brown and 
Brennand [19], but show the following behavior: 
(a) At any temperature the hysteresis loss increases 
with induction, following the Steinmetz law W=aB n 
quite well; (b) at any induction the hysteresis loss 
increases with decreasing temperature; (c) the rela- 
tive increase in the loss is independent of induction: 
at 20° K, for instance, the loss appears to have in- 
creased by 20 percent over the loss at ambient 
temperature for all inductions between 2 and 10 
kilogauss. 



Among those contributing to the work of this study 
were: W. F. Mayer (guest worker from AiResearcli 
Manufacturing Co., summer 1956) to whom credit is 
due both for the basic design of the turbine for the 
bearing tester, and guidance in the initial stages of 
the study; R. J. Corbett (guest worker from Beech 
Aircraft Corp., summer-fall 1956) who assisted in 
the detailed design of the bearing tester and was also 
very helpful in comments and criticism of this 
report; the following part-time student employees: 
D. W. Murphy (who carried out most of the sonic 
velocity calculations), J. R. Brennand (who per- 
formed most of the hysteresis loss tests), and a 
number of others who assisted on occasion. 
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